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We analyze the properties of tlie contact process witii long-range interactions by the use of a 
kinetic ensemble in which the total number of particles is strictly conserved. In this ensemble, both 
annihilation and creation processes are replaced by an unique process in which a particle of the 
system chosen at random leaves its place and jumps to an active site. The present approach is 
particularly useful for determining the transition point and the nature of the transition, whether 
continuous or discontinuous, by evaluating the fractal dimension of the cluster at the emergence of 
the phase transition. We also present another criterion appropriate to identify the phase transition 
that consists of studying the system in the supercritical regime, where the presence of a "loop" char- 
acterizes the first-order transition. All results obtained by the present approach are in full agreement 
with those obtained by using the constant rate ensemble, supporting that, in the thermodynamic 
limit the results from distinct ensembles are equivalent. 
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I. INTRODUCTION 

The study of equivalence of ensembles that describe 
systems in thermodynamic equilibrium, the Gibbs en- 
sembles, has a long tradition in statistical mechanics. 
Thermodynamic properties evaluated in distinct ensem- 
bles will be the same if they are equivalent. In the case of 
equilibrium systems, the ensembles are set up according 
to well known prescriptions that take account the exis- 
tence of a Hamiltonian P, i]. The possibility of repre- 
senting nonequilibrium systems in distinct ensembles was 
first considered by ZifF and Brosillow [sj when they used 
the constant coverage ensemble to study first-order tran- 
sitions in a catalytic reaction model [j , originally defined 
in the ensemble of constant rate. Afterwards, Tome and 
de Oliveira 0] introduced the conserved contact process, 
a version of the contact process 0] in the ensemble 
of constant particle number. The equivalence between 
the constant rate and constant particle number ensem- 
bles for the contact process in the stationary regime has 
been proved by Hilhorst and van Wijland in the thermo- 
dynamic limit [8]. Later on, the equivalence of ensembles 
in the stationary regime was extended for other systems 
with short-range interactions with distinct annihilation 
dynamics (9l.[To|. 

In equilibrium statistical mechanics, the equivalence 
of ensembles is granted for homogeneous and extensive 
systems in the thermodynamic limit P, 0, [lH ■ On the 
other hand, for non-extensive systems, the canonical and 
microcanonical ensembles are not equivalent . Exam- 
ples of physical problems in which canonical and micro- 
canonical ensembles are not equivalent are nuclei, atomic 
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clusters in which the range of forces between their con- 
stituents are comparable to the system size, systems at 
the first-order transition with phase separation [Til . [Tsj 
and models with long-range forces [IlQIil- In partic- 
ular, this latter set of problems has deserved a great inter- 
est not only in statistical mechanics, but also in other ar- 
eas, such as nuclear physics, astrophysics [l^ and plasma 
physics . Nonequilibrium systems with long-range in- 
teractions have also been proposed and in this case, we 
may ask whether the equivalence of ensembles may be 
valid for such systems. 

In the context of nonequilibrium statistical mechanics, 
systems displaying long-range interactions have been pro- 
posed originally as more realistic models describing the 
spreading processes, instead of short-range systems. In 
particular, an anomalous model of the directed percola- 
tion has been proposed by MoUison In this problem, 
the infection probability performs a Levy flight decaying 
with the distance r as a power- law relation l/r'^+'^, where 
d is the spatial dimension of the system and a is control 
parameter. The critical behavior of anomalous directed 
percolation displays a family of universality classes that 
have been studied by Grassberger ^18|, by Janssen et al 
[l^ who have considered field-theoretic renormalization 
group, by Hinrichsen and Howard [20| who performed 
numerical simulations and recently by Tessone et al [2lj 
who analyzed a class of spatially extended chaotic sys- 
tems with power-law decaying interactions. In all cases 
above, the critical exponents vary continuously with the 
parameter a. More recently, another class of nonequilib- 
rium systems with long-range interactions, named con- 
tact processes with long-range interactions, have been 
introduced by Ginelli et al. [13, [H, [23|, inspired by 
pinning-dcpinning transitions in nonequilibrium wetting 
phenomena. By varying the control parameters, the con- 
tact process with long-range interactions exhibits a richer 
phase diagram than the usual short-range contact pro- 
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cess, with distinct universality classes [23| and discontin- 
uous phase transitions p^ . [23j . 

Two aims concern us in this paper. First, we wish 
to analyze the equivalence of ensembles in nonequilib- 
rium systems with long-range interactions. The equiv- 
alence of ensembles has been studied in systems with 
short-range interactions [1, [l^, [2^, H^, but not in sys- 
tems with long-range interactions. We will consider here 
the (T— contact process with long-range interactions, in- 
troduced by Ginelli et al [12] ■ Second, we wish to show 
that the present approach offers a new procedure for an- 
alyzing first-order transitions of systems with absorbing 
states. Due to the existence of the absorbing states, stan- 
dard procedures used successfully in the study of dis- 
continuous equilibrium phase transitions may not work 
well when applied to describe nonequilibrium first-order 
transitions. In the present case, the first evidence for a 
discontinuous transition will be given by the presence of 
"loops" in the rate versus density curve, as long as the 
system is finite. As it will be shown, the loops disappear 
in the thermodynamic limit giving rise to a tie line, the 
true signature of a first-order transition. Another ad- 
vantage of using the present approach is verified when 
one analyzes the system at the emergence of the phase 
transition, since the nature of the phase transition can be 
inferred simply by examinating the structure of particles. 
As we shall see, if, at the emergence of the transition, the 
cluster is fractal, the transition will be continuous; if the 
cluster is compact, the transition will be discontinuous. 
A compact cluster provides us an evidence of a discon- 
tinuous transition. 



II. MODEL 

A. Constant rate ensemble 

The one dimensional cr— contact process with long- 
range interactions [2^ is defined in a chain of L sites with 
periodic boundary conditions as follows. To each site i of 
a one-dimensional lattice is attached an occupation vari- 
able iji that takes the values or 1 according whether 
the site i is empty or occupied by a particle, respectively. 
The process is composed of spontaneous annihilation of 
a single particle (1^0) and catalytic creation of a par- 
ticle (0^1). Particles are created only in empty sites 
in which at least one of its nearest neighbor sites is occu- 
pied by a particle. These empty sites are named active 
sites. The rate of creation depends on the length of the 
island of empty sites next to the active site and decreases 
algebraically with the size of the island, introducing an 
effective long-range interaction. The total transition rate 
'Wi{r]) is given by 

+ (1) 

where a is a parameter and 

^nv)=V^, (2) 



describes a spontaneous annihilation and 
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+ 2 + j^)Vi+iVtVi-i ■ ■ ■ Vi-t+m~e, (3) 

i=i 

describes a catalytic creation that depends on the length 
£ of the island of empty sites, where a and a are pa- 
rameters and we have considered the shorthand notation 
fji = 1 — r]i in Eq. ([3]). When a = one recovers the 
original short-range contact process @, 0| ■ 

Here, numerical simulations in the constant rate en- 
semble is performed as follows. A particle is chosen at 
random from a list of occupied sites. It is annihilated 
with probability p — a/ (l+a+a). With probability 1—p, 
a new particle may be created next to the chosen parti- 
cle. This is done by choosing first one of its two nearest 
neighbors with equal probability and then a particle will 
be actually created with probabihty q = {l+ai~'^) / {1+a) 
provided the chosen site is empty. This algorithm gives 
the following ratio between the creation and annihilation 
of particles {l—p)q/2p = (l + a£~'^)/2a, which is equiva- 
lent to that considered by Ginelli et al [22] , in which the 
creation and annihilation of particles occur with rates 
A(l -I- a£~"') and 1, respectively as long as a is related to 
A by a — 1/2A. The increment of the time is given by 
1/Np, where Np is the number of occupied sites. 

For large values of a, the system is constrained into the 
absorbing state, in which no particles are allowed to be 
created. Decreasing the parameter a, a phase transition 
to an active state takes place, whose location depends on 
the parameters a and a. For a fixed value of a (we take 
here to be a = 2, as considered by Ginelli et al. [2^) and 
a > 1 the transition is continuous with critical exponents 
belonging to direct percolation (DP) universality class. 
For < a < I the transition becomes discontinuous. 
Ginelli et al found that, the crossover between the 
discontinuous and second-order phase transitions occurs 
at cr = 1.0. 



B. Constant particle number ensemble 

In the constant particle number ensemble, the control 
parameter is the total particle number n. Again we used 
a chain with L sites with periodic boundary conditions. 
Particles are neither created nor annihilated. Instead, a 
particle leaves their place and jump to an empty site. 
However, the jump process is not unrestrictive process, 
since their occurrence must be consistent with the rules of 
their version in the constant rate ensemble. More specif- 
ically, the dynamics that characterizes this ensemble is 
defined as follows. A particle of the system, chosen at 
random, leaves its place, located for example at the site 
i, and jumps to an empty site placed at the site j, also 
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chosen at random. The jump rate will depend on the 
specific rule of the considered model. One may define 
the jump process by a dynamics in which both creation 
and annihilation occurs simultaneously, according to the 
transition rate (77) is given by U,um 



(4) 



To see that this transition rate leads to a dynamic that 
is equivalent to that given by Eq. |T]) let us consider 
the total rate ^iWij{r])/L = J2i^ii''l)^j{v)/L in which 
particles jump to site j. In the thermodynamic limit 
the system size L ^ 00 and the sum X)i '-^fl'?)/-^ ^P" 
proaches {uj^{rj)), by the law of large numbers, so that 
^iWij{ri)/L = {ijjf{r]))ujj{r]). By an analogous argu- 
ment the total rate in which particles leave the site i is 
J^j^ijiv)/!^ = {ui'j{r]))u;f{r]). Comparing with the Eq. 
([T]), the averages {ujf{r])) and (cuKr])) should be propor- 
tional to 1 and a, respectively, thus 



K(^)) 



(5) 



where ujj{r]) and ^"(77) are given by Eqs. ^ and ([3]), 
respectively, so that Eq. ^ is given by 



P 



(6) 



where {uj°i{rj)) = {rji) = pis the density of particles. This 
formula allow us to evaluate the parameter a within re- 
spect to the ensemble of constant particle number. The 
average {u!j{r])) can be understood as "density of active 
sites" , that is, the density of empty sites in which parti- 
cles may jump to them. It can be evaluated directly from 
numerical simulations by computing Eq. ([3]) that is zero 
and (1 -I- ai~'^)/2 according whether the site i is occu- 
pied or empty, respectively, where £ is the length of the 
island of inactive sites in which the active site is located. 
For time-dependent regime, however, Eq. ^ might not 
always hold. For instance, if the initial state is such that 
averages {uJj{7])) and (w°(ry)) are not constants, Eq. ^ 
cannot be satisfied [1, Q . 

The actual numerical simulation of the constant par- 
ticle number ensemble is realized as follows. An empty 
site surrounded by at least one particle (active site) is 
chosen at random. Next, a particle of the system, also 
chosen at random, jumps to the active site with probabil- 
ity {l + a£^'^)/{l + a). The constant factor l/(l-|-a) 
is used in order to guarantee that pi < 1, since 1 -|- a£~^" 
may be greater than 1. In the constant rate ensemble, 
the control parameter may be included in the probabil- 
ity of choosing each subprocess (creation or annihilation). 
However, in the present case, as both processes occur si- 
multaneously and the parameter a is not constant, this 
procedure is not possible. 

In contrast to the constant rate ensemble, the con- 
served contact process does not have, strictly speaking, 
an absorbing state. This fact constitutes an useful tool in 



the study of phase transitions because there is no danger 
of falling into the absorbing state as happens in numer- 
ical simulations of the constant rate ensemble. In the 
conserved ensemble (constant particle number ensemble) 
the equivalent of the absorbing state is the state with 
zero density, as is the case of an infinite system with a 
finite number of particles, also named subcritical regime. 



III. CLUSTER APPROXIMATIONS 

The first analysis we have performed here consists of 
studying the cr— contact process by means of cluster ap- 
proximations. We will consider here approximations at 
the level of one and two sites. 

At the level of one-site approximation, we use the fol- 
lowing approximation for the probability of a string of 
sites: 



(7) 



In this case, it suffices to consider the dynamic variables 
P{1) = p and -P(O) — 1 ~ p. In the steady state, one has 
the following relation 



a = 1 — p + ap 



[i-pY 



(8) 



e=i 



This equation has already been obtained by Ginelli et al 
[13, m. For < (7 < 1, Eq. ^ establishes a discontinu- 
ous transition between an absorbing and an active state 
that can be viewed by the existence of "loops" [23| • As it 
will be seen later, numerical simulations in the constant 
particle number ensemble also presents "loops" for finite 
systems in the interval < a < 1. However, Eq. ([S]) 
also exhibits loops for cr > 1, in contrast to numerical 
simulations. 

In order to obtain improved results, we have derived 
relations by considering correlations of two sites. In this 
case, the probability of a string of sites is approximated 
by 



P{rii,V2)PiV2,m)---P{rii-i,m) 
P{m)P{V3)-P{ve-i) 



(9) 

The dynamic variables are now P(1),P(0) and nearest- 
neighbor joint probabilities P(11),P(10) ~ P(01) and 
P(00). Only two of them are independent, so that it is 
necessary to write down two equations. By solving these 
equations numerically, we find the behavior of a versus 
p for several values of cr, as showed in Fig[T] In contrast 
to the previous case, two-site approximations supports a 
change in the nature of the phase transition by increasing 
the parameter cr, passing from first-order to second-order 
at CT = (1.3 ± 0.1). This can be identified by disappear- 
ance of loops. Although the results obtained from two- 
site approximations agree qualitatively with numerical 
simulations, the location of the transition points as well 
as the value of a that characterizes the crossover between 
first-order and second-order transitions are incorrect. 
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FIG. 1: Quantity a versus p for several values of a at the level 
of two-site mean field approximation. A crossover from a first- 
order to a continuous transitions occurs ai a — (1.3 ± 0.1). 
Symbols have been used in order to distinguish the distinct 
values of a. 



TABLE I: Results of numerical simulations coming from 
the constant particle number ensemble (second and third 
columns) and from the constant rate ensemble (fourth and 
fifth columns) for a system size L = 20000. 



2.0 0.7500 0.2964(4) 0.29638 0.7509(9) 

1.5 0.6900 0.33481(1) 0.33480 0.6900(1) 

1.0 0.7000 0.37007(1) 0.37000 0.7004(2) 

0.5 0.8500 0.32887(8) 0.32890 0.85008(8) 

0.4 0.9000 0.24948(2) 0.24948 0.90000(1) 
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IV. NUMERICAL RESULTS 



Numerical simulations of the cr— contact process were 
performed for a = 2 and several values of a. We have 
used 2 X 10'' Monte Carlo steps to evaluate the appro- 
priate quantities, after discarding a sufficient number of 
steps to reach the stationary state. Here one Monte 
Carlo step corresponds to n jumping processes. In or- 
der to compare results obtained from distinct ensembles 
we have simulated the contact process with long-range 
interactions not only in the constant particle number en- 
semble, but also in the constant rate ensemble. We have 
first determined the quantity a in the conserved (con- 
stant particle number) ensemble, by using Eq. ([T]), for 
several densities p. Next we used the values of a = a 
to perform numerical simulations in the constant rate 
ensemble which in turn gives us the average density p. 
According to Table I, the excellent agreement between 
the results confirms the equivalence of ensembles. How- 
ever, numerical simulations provide distinct results at the 
phase coexistence, even when one considers large system 
sizes. For example, for L = 20000 and a = 0.5, simula- 
tions in the constant rate ensemble for a — 0.42313 gives 
p = 0.76600(4). This value of a corresponds to two den- 
sities in the constant particle number ensemble p — 0.60 
and p = 0.766. As it will be shown later, both ensem- 
bles become equivalent at the phase coexistence in the 
thermodynamic limit. 

From now on we will drop the bars over a and p. 



FIG. 2: Quantity a versus 1/n for the subcritical regime. 
The left, center and right panels correspond to a = 2, a — 0.5 
and a — 0.4, respectively. The straight lines fitted to the 
data points give the extrapolated critical points ao = 0.35497, 
ao = 0.41389 and ao = 0.41714. 

A. Subcritical regime 

Let us consider a finite number of particles placed on 
infinite lattice. In this situation the density of particles 
vanishes and the system is constrained to remain in the 
subcritical regime. The actual simulation of the subcrit- 
ical regime is done by using a finite lattice and check 
whether a particle reaches the border. If a particle never 
reaches the border, the system size may be taken as infi- 
nite. The size of the system was chosen to be big enough, 
so that no particles have reached the border for the max- 
imum time considered. An important feature of the sub- 
critical regime is that the increase of the total number of 
particles makes the system approach the transition point 
at a = OfQ. The expected value of a obtained for a fixed 
number of particles n a ppr oaches its asymptotic value ao 
according to [1, ^ ^JQ^ 

a — ao ^ —. (10) 
n 

A linear extrapolation of a versus 1/n when n ^ oo 
gives ao- In Fig- [H we have plotted the numerical val- 
ues of a determined from simulations of an infinite sys- 
tem with n particles for different values of a. Numer- 
ical extrapolations obtained by using Eq. (jTU]) give us 
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FIG. 3; Typical run for the conserved long-range contact pro- 
cess with n = 50 particles for a = 2 starting from a random 
configuration. A unit of time, or a Monte Carlo step, corre- 
sponds to n particle jumps. We have discarded 3 x 10^ initial 
Monte Carlo steps. 



FIG. 4: Typical run for the conserved long-range contact pro- 
cess with n = 50 particles for o — 0.5 starting from a random 
configuration. A unit of time, or a Monte Carlo step, corre- 
sponds to n particle jumps. We have discarded 3 x 10^ initial 
Monte Carlo steps. 



ao = 0.41714(6), ao = 0.41389(1) and ao = 0.35497(6) 
for a = 0.4, a = 0.5 and a = 2, respectively, which agree 
very well with estimates ao = 0.4172(1), ao = 0.41382(3) 
and ao = 0.3548(2), obtained from the constant rate en- 
semble. The same agreement is verified for other values 
of CT, revealing the utility of the present procedure for lo- 
cating the transition point. In the section IV C, we will 
show results from time-dependent numerical simulations 
for a = 0.4. For a = 0, the phase transition takes place at 
ao = 0.303228(2V as expected for the usual short-range 
contact process @. 

The first procedure considered here for classifying the 
phase transition consists of examinating the spatial struc- 
ture of particles at the emergence of the phase transition. 
In Figs. [3] and [3] we have plotted the time evolution for 
two distinct values of a. For o- = 2, the system generates 
fractal clusters, whereas for a = 0.5 they are compact. 

A measure of the size of the cluster of particles is given 
by the quantity R given by maximum distance between 
two particles of the cluster [2^, [2^. As long as n is finite, 
R is also finite, but it diverges when n — > oo. It is related 
to the total particle number n through the relation 

Rr^n^/"^^, (11) 

where dp is the fractal dimension [2^ [2^. The quan- 
tity R will be evaluated here by measuring the end-to- 
end spread of the cluster [13, [M US HI- In Fig- 4 
we show the log-log plot of R versus n in the subcriti- 
cal regime for several values of a. For cr > 1, the slopes 
of the straight lines fitted to the data points are consis- 
tent with 1.33704(4), as expected for the DP universal- 
ity class whose clusters generated at the criticality [2^ 
have fractal dimension dp = 0.74792(2) [fj], whereas for 
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FIG. 5: Log-log plot of the average cluster size R versus the 
total number of particles n for several values of a at the sub- 
critical regime. The straight lines from top to bottom have 
slopes 1.33704 and 1, respectively. The data points have been 
shifted in order to avoid overlapping. 



< cr < 1 the slopes are equal to 1. The change in the 
fractal dimension shows the change in the nature of the 
phase transition, passing from continuous to discontinu- 
ous, as one decreases a. A fractal dimension equal to the 
Euclidean dimension is actually the signature of a first 
order transition, since we have a compact cluster that 
remains finite in the thermodynamic limit. To see this, 
let us evaluate the density p* of the cluster of particles 
p* — n/ R. It is worth emphasizing that although the to- 
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tal density p = n/Loi the system is zero, since the system 
is constrained in the subcritical regime, the quantity p* 
is finite. If the phase transition is first order, p* will have 
a non-zero value in the thermodynamic limit. Rewriting 
Eq. pT|) in terms of p*, one has the following relation 
p* = n/R - Therefore, for < ct < 1, 

where dp = 1, p* — > Po' when n — > cx), corresponding 
to the active phase pg in coexistence with the absorbing 
phase (p = 0). The density Pq is determined simply by 
the inverse of the slope of curves in Eq. (fTTj) . In contrast, 
for a second order transition in which dp < 1, p* vanishes 
when n ^ oo, because p* — n/R ^ n-(i-dF)/dF _^ q 
when n — > cxD. The change in the nature of the phase 
transition for a close to 1 is in agreement with results 
by Ginelli et al [22*1 who used time-dependent numerical 
simulations and distribution of inactive islands to iden- 
tify the order of the phase transition. 

We close this section by comparing the density pg of 
the compact cluster obtained from both ensembles. As 
it was mentioned above, in the constant particle num- 
ber ensemble the density pg is evaluated directly, since 
it is simply the inverse of the slope of Eq. (fTTj) . To 
determine pg in the constant rate ensemble, we use the 
procedure adopted by Dickman [s^, HH that consists of 
dividing the system into blocks of 100 sites and deter- 
mining histograms of such blocks density profiles at the 
phase coexistence. For a — 0.5 at the phase coexistence, 
the histogram show a bimodal distribution with a peak 
at p* = and another one at pj = 0.77. This agrees 
with the density of the compact cluster pg — 0.777(1), 
obtained from the present approach (constant particle 
number ensemble). For other values of a in the inter- 
val < a < 1, the densities of compact clusters at the 
phase coexistence obtained from both ensembles agree 
very well. 



B. Supercritical regime 

The supercritical regime is characterized by nonzero 
values of the density p. To simulate the cr— contact pro- 
cess in the supercritical regime, we have considered lat- 
tice sizes ranging from L = 200 to i = 20000. In Figs. [5] 
and [7] we have plotted the quantity a versus p for several 
values of L for a = 2 and a = 0.5, respectively. The 
curves exhibit a strong dependence on L that is qualita- 
tively different whether we consider (T>lor0<cr<l. 
In the first case (exemplified in Fig. [S] by u — 2), the 
curves are strictly decreasing functions and cumulated, 
when L — > oo, into a strictly decreasing function. 

In the second case (exemplified in Fig. [7] by cr = 0.5), 
the curves are no longer decreasing and present "loops" , 
that can be associated to a coexistence of two phases. 
As it was shown previously, the occurrence of "loops" in 
numerical simulations agrees qualitatively mean-field re- 
sults. However, in contrast to the mean-field approach, 
the existence of "loops" in numerical simulations is a fi- 
nite size effect that disappears when L —>■ oo, giving rise 
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FIG. 6: Quantity a versus p for several values of the system 
size L for a — 2. The inset shows the collapse of the data 
by using relation ((12} and (|13|) where y = L^^^^la — ao\ and 
X = L^^'^^ p for the supercritical regime and y = L''^ |a — ao| 
and X = L''~'*^p for the subcritical regime. 
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FIG. 7: Quantity a versus p for several values of the system 
size L for cr = 0.5. The horizontal straight line is the tie 
line at Q = Qo = 0.41389. The dashed vertical line indicates 
the density p* = 0.777 of the coexisting active phase. The 
inset shows the collapse of the data by using relation (|14p 
y = (q - ao)L'' . 

to a horizontal tie line that connects the two coexisting 
phases. 

For CT > 1, where the phase transition is second order, 
the finite size scaling for the density is given by Q 

p = L-'3/--/(eLi/^-), (12) 

valid for the supercritical regime, where e = a — ao and (3 
and are critical exponents associated to the order pa- 
rameter and the spatial correlation length, respectively; 
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FIG. 8: Log-log plot of a — qq versus L for several densities p 
for G = 0.5 (left panel), g = 0.4 (center panel) and g = 0.25 
(right panel). From left to right, the straight lines at each 
panel have slopes 0.5, 0.4 and 0.25. The data points have 
been shifted in order to avoid overlapping. 

and H 

p = L-'^+''^f{eL'^''), (13) 

valid in the subcritical regime, where f{x) is an universal 
function. A data collapse obtained by using the finite size 
scalings (|12p and and the best estimates of the DP 
critical exponents (3 = 0.276486(8), = 1.096854(4) 
and dp = 0.74792(2) % is shown in the inset of Fig. M 

On the other hand, for < u < 1, where the phase 
transition is first order, Eqs. (fT^ and (fT5|) are not valid. 
To relate the parameter a with the system size L, we 
assume the following behavior of a in the region of the 
phase coexistence 

a-ao^L-", (14) 

where ao is the value of a at the coexistence, or the value 
of a where the tie line is located. 

An argumentation concerning Eq. (|14|) is given by as- 
suming the one-site mean-field approximation given by 
Eq. (IH). For < cr < 1, the sum on the right side of 
Eq. ^ can be replaced by an integral, which to leading 
order in p, reduces to [2^. l23j 

a — ao — —p + aT{l — (j)p'^ . (15) 

where ao = 1 is its one-site mean-field transition value 
and T{x) is the gamma function. Relating p with the 
system size L, we have that a — L~'' ^ because the 
second term in the right-hand side of Eq. (|15p dominates 
over the first. 

The correctness of Eq. ([14]) can be checked by the 
log-log plot of a — ao versus L by using the estimation 
of ao available from the study of the subcritical regime. 
Indeed, as shown in Fig. [5] the slopes of the straight 



FIG. 9: Dependence of the parameter a on L~'' for g = 0.5 
(left) and g = 0.4 (right) considering several densities p. The 
full circles are extrapolated values of qq obtained from the 
subcritical regime. 



lines fitted to the data points for several densities are 
consistent with the values of a. We have repeated the 
analysis for other values of a in the interval < cr < 1, 
whose dependence of a on the system size L at the phase 
coexistence is also described by Eq. ((HI). For higher 
densities close to p* , in which one observes a peak of a 
versus p, it is necessary to consider larger system sizes in 
order to reach the asymptotic behavior described by Eq. 

as can be seen for the highest densities in the left 
and center panels of Fig. [51 respectively. 

A complementar analysis, but fully equivalent, con- 
sists in assuming Eq. (|14p and using it for determining 
estimates of ao by numerical extrapolation, as shown in 
Fig. O The agreement between extrapolated ao using 
Eq. ITU) and those estimates of ao available from the 
subcritical regime and spreading simulations (reported in 
the next section) confirms the equivalence of ensembles 
in the thermodynamic limit at the phase coexistence. 

As a final checking, we plot in the inset of Fig. [7] a col- 
lapse of the data for different system sizes by considering 
the variable y = {a — OLf))L" versus p. This confirms once 
more the conjecture given by Eq. (|14p . 



C. Time dependent numerical simulations 

Here, we show explicitly results from time-dependent 
numerical simulations for the cr— contact process consid- 
ering cr ~ 0.4, in order to compare the results for the 
estimation of ao from distinct ensembles. 

Starting from a configuration close to the absorbing 
state, this procedure consists in determining the time 
evolution of appropriate quantities, such as the survival 
probability Ps{t), the total number of particles Np{t) and 
the mean square spreading (t) of the active region. At 
the emergence of a second-order transition, these quan- 
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t t 

FIG. 10: Temporal evolution of the total number of particles 
Np (right) and the survival probability Ps (left) for the ordi- 
nary a— contact process for a = 0.4. From top to bottom, we 
have Q = 0.4167, 0.4171, 0.4172 and 0.4176. 

titles are described by following power-law behaviors 

Ps{t)--t-^, Np{t)^f' and R^{t)r^e, (16) 

where 5, rj and z are their associated critical exponents. 
They are related to the fractal dimension introduced in 
Eq. PT|) through the relation dp = 2{r] + 6)/z. 

Although the order parameter presents a jump in a 
nonequilibrium first-order transition, some dynamic vari- 
ables are also characterized by dynamic exponents [s^ . 
In the present case, the quantities from Eq. p6p are 
expected to be described by the same exponents as the 
Glauber-Ising model at zero temperature, which values 
are S = 1/2, 77 = and z = 1 [H]. Relating these 
exponents, we have that dp = 1, in agreement with re- 
sults obtained previously from analysis in the subcritical 
regime for < cr < 1. 

In Fig. ITUl) . we show the plot of temporal evolu- 
tion of the quantities Pg, Np for some values of a. For 
a = 0.4172, the quantities Np{t) and Ps{t) follow indeed 
a power-law behavior described by Eq. (|16p . whose expo- 
nents are consistent to rj = and S = 1/2, respectively. 
This estimation of ao agrees very well to those obtained 
from the constant particle number ensemble, confirming 
the equivalence of ensembles. 

V. CONCLUSION 

In this paper we have studied a nonequilibrium model 
with long-range interactions, named cr— contact process, 
in the ensemble of constant particle number. The equiv- 



alence of ensembles is confirmed by the excellent agree- 
ment between the numerical results coming from both 
ensembles. All results obtained by the present approach 
are in full agreement with those obtained by Ginelli et al 
[2^ . We believe that the present approach may be partic- 
ularly useful for studying discontinuous phase transition, 
since it is possible to identify the nature of the transi- 
tion by measuring the spatial structure of particles at 
the transition. The interesting feature of the present ap- 
proach concerns the study of an infinite system with a 
finite number of particles. As the number of particles 
is increased, the systems naturally approaches the criti- 
cal point and, in this sense, it behaves like self-organized 
critical systems. 

Another advantage of the present approach is verified 
when one studies the system in the supercritical regime. 
The dependence on the system size is distinct in both 
regimes as shown in Figs. [S] and [T] The presence of 
a "loop" is a strong indication of a first-order transition. 
Of course, results coming from finite systems, such as the 
"loop" presented in Fig. [7] is a particularity of the con- 
served ensemble, that disappears in the thermodynamic 
limit. 

We remark that the use of this procedure in the con- 
stant rate ensemble is not possible since numerical simu- 
lations of finite systems will present a jump in the order 
parameter close to the transition point, even in the case 
of a second-order transition, due to the presence of the 
absorbing state. Recently, de Oliveira and Dickman [s^l 
have proposed a method, named quasi-stationary simula- 
tions, that improves the accuracy of results, even at the 
emergence of the phase transition. It consists in simu- 
lating the system in the constant rate ensemble in the 
standard way. However, whenever the system enters in 
the absorbing state, a non-absorbing configuration is cho- 
sen from a list of saved periodically configurations and 
then system returns to the active state. This method 
has revealed an useful tool in the study of systems with 
absorbing states [1^ [3l|, [33|. Other techniques, such 
as hysteretic analysis in a constant coverage ensemble 
have also been used to study discontinuous transition in 
nonequilibrium systems 0, [H, [11]. These approaches 
are, however, distinct from ours, in the sense that the 
standard ensemble is initially used to generate a station- 
ary configuration and then the system is switched to the 
constant coverage ensemble. 
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